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Q , ABSTRACT 

^ ■ Since the 1980's it has been becoming increasingly clear that the Solar System's 

irregular satellites are collisionally evolved. The current populations are rem- 
04 . nants of much more massive swarms that have been grinding away for billions 

<N ■ of years. Here, we derive a general model for the collisional evolution of an ir- 

regular satellite swarm and apply it to the Solar System and extrasolar planets. 
p ^ ' . The model uses a particle in a box formalism and considers implications for the 

size distribution, which allows a connection between irregular satellite popula- 
tions and predicted levels in the resulting dust cloud. Our model reproduces the 

■ Solar System's complement of observed irregulars well, and suggests that the 

competition between grain-grain collisions and Poynting-Robertson (PR) drag 

q , helps set the fate of the dust. In collision dominated swarms most dust is lost 

?-h ■ to interplanetary space or impacts the host planet, while PR dominated grains 

c/j | spiral in towards the planet through the domain of regular satellites. Because 

Ci . swarm collision rates decrease over time the main dust sink can change with 

time, and may help unravel the accretion history of synchronously rotating reg- 
ular satellites that show brightness asymmetries, such as Callisto and Iapetus. 

■ Some level of dust must be present on AU scales around the Solar System's giant 
(X) ' planets if the irregular satellites are still grinding down, which we predict may 
I/"") , be at detectable levels. We also use our model to predict whether dust produced 
OO ■ by extrasolar circumplanetary swarms can be detected. Though designed with 

planets in mind, the coronagraphic instruments on JWST will have the ability 
to detect the dust generated by these swarms, which are most detectable around 
planets that orbit at many tens of AU from the youngest stars. Because the col- 
lisional decay of swarms is relatively insensitive to planet mass, swarms can be 
much brighter than their host planets and allow discovery of Neptune-mass plan- 
ets that would otherwise remain invisible. This dust could have been detected 
by HST ACS coronagraphic observations, and in one case dust may have already 

■ been detected. The observations of the planet Fomalhaut b can be explained as 
j scattered light from dust produced by the collisional decay of an irregular satel- 
lite swarm around a ^10 planet. Such a swarm comprises about 5 Lunar 
masses worth of irregular satellites. Finally, we briefly consider what happens if 
Fomalhaut b passes through Fomalhaut's main debris ring on a coplanar orbit, 
which allows the circumplanetary swarm to be replenished through collisions 
with ring planetesimals. This scenario, in which the planet is at least of order an 
Earth mass, may be ruled out by the narrow structure of the debris ring. 

Key words: circumstellar matter - stars: planetary systems: general - Solar 
System: planets and satellites 



1 INTRODUCTION 

With a penchant for retrograde, barely-bound, high- 
Email: gkennedy@ast.cam.ac.uk eccentricity orbits and flatter-than-usual size distribu- 
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tions, the irregular satellites are one of the Solar Sys- 
tem's rebel populations. Their presumed capture into 
these unusual orbits around the Solar System's giant 
planets has long been a puzzle. Because irregulars ex- 
ist at all four outer planets, capture mechanisms specific 
to formation of gas giants like Jupiter and Saturn (e.g. 
Heppenheimer & Porco 1977; Pollack et al. 1979) are not 
general enough. Dynamical mechanisms, which do not 
rely on the existence or growth of large gaseous atmo- 
spheres, can be applied to gas and ice-giants alike and 
are therefore favoured (e.g. Colombo & Franklin 1971). 
Nesvorny et al. (2007) recently proposed a different dy- 
namical mechanism as part of a unified model of outer So- 
lar System formation (Gomes et al. 2005; Tsiganis et al. 
2005), where irregulars are captured during a period of 
instability via planet-planet interactions. 

An apparent weakness of the planet-planet inter- 
action mechanism is that the irregular satellites should 
share the same size distribution as Jupiter's Trojan as- 
teroids, which in this model were captured at the same 
time from the same population (Morbidelli et al. 2005). 
In fact the differences in the Trojan and irregular satel- 
lite size distributions are marked, with irregular satellites 
being much flatter for sizes larger than about 10 km. To 
overcome this hurdle, Bottke et al. (2010) showed that 
while the size distribution was indeed initially steeper 
like the Trojans, 4.5 Gyr or so of collisional evolution is 
sufficient to reduce primordial irregular satellite popula- 
tions to a size distribution that matches those currently 
observed. This result implies that there were previously 
more irregular satellites and perhaps most significantly, 
that copious amounts of dust were produced during the 
depletion of these satellites. 

The evidence for collisional evolution of irregular 
satellites has been mounting for some time. Kessler 
(1981) showed that the four prograde irregulars known 
to orbit Jupiter at the time had a relatively short col- 
lisional lifetime. The advent of large-format CCD sur- 
veys since the turn of the century has seen a dra- 
matic increase the number of irregulars and made fur- 
ther theoretical advances possible (e.g. Gladman et al. 
2001; Sheppard & Jewitt 2003; Holman et al. 2004; 
Jewitt & Haghighipour 2007). Nesvorny et al. (2003) 
noted that irregular satellites around planets closer to the 
Sun have larger orbits (in Hill radii) , suggesting that the 
lack of satellites closer to Jupiter is due to their erosion 
through collisions, which proceed at a faster rate closer to 
the planet. Using numerical integrations to derive average 
orbital elements, they also proved the existence of colli- 
sional families (see also Gladman et al. 2001). This latter 
discovery is particularly important, because the current 
collision rate amongst the irregular satellites is in some 
cases too low to explain their existence (Nesvorny et al. 
2003). The inference is again that the number of irregular 
satellites, and thus their collision rate, was much greater 
in the past and that they decayed through collisions to 
the current level. 

Coinciding with these theoretical advances was the 
discovery by direct imaging of Fomalhaut b (Kalas et al. 
2008), an extrasolar planet predicted to exist based on 
the elliptical orbit of Fomalhaut's circumstellar debris 
ring (Kalas et al. 2005; Quillcn 2006). The ring structure 



suggests that the planet is less than 3Mj up , though the 
planet could be much less massive (Chiang et al. 2009). 
This discovery appears unrelated to irregular satellites, 
but the inability of Hubble, Keck, and Gemini photome- 
try to pin down whether the planet looks like a planetary 
atmosphere or reflected starlight provides the link. 

While the planets discovered to orbit HR 8799 ap- 
pear to be consistent with ~1000K substellar-mass ob- 
jects (Marois et al. 2008), multi- wavelength photome- 
try of Fomalhaut b appears bluer than expected for a 
200 Myr old gas-giant planet. Specifically, Fomalhaut b 
has so far defied detection at wavelengths longer than 
1/mi, leading Kalas et al. (2008) to suggest that the spec- 
trum is actually starlight scattered from an optically 
thick circumplanetary disk of about 20 Jupiter radii. 
Though such a scenario is plausible, we argue that dust 
produced by a swarm of colliding irregular satellites is 
also a possibility. 

Given that the Solar System's irregular satellite com- 
plement decayed to its current state through collisions 
and the exciting possibility that Fomalhaut b may har- 
bour the first circumplanetary dust seen outside the So- 
lar System, the time seems right to consider whether such 
clouds of irregular satellites could be visible around extra- 
solar planets. In the following sections, we derive a sim- 
ple model for the evolution of a circumplanetary satellite 
swarm and the all-important extrasolar observable — the 
dust. We compare our model with the Solar System irreg- 
ulars, and comment on the fate and observability of dust. 
We then apply the model to circumplanetary swarms 
around extrasolar planets. Finally, we explore what kind 
of satellite swarm could exist around Fomalhaut b, and 
what constraints this proposed swarm puts on the planet 
mass. 



2 MODEL OF A CIRCUMPLANETARY 
SWARM 

The irregular satellite swarms described in this paper 
have not knowingly been detected around other planets 
that orbit other stars. Therefore, like the pre-1995 days 
of planet formation theory, we must take cues from the 
Solar System. However, based on experience gained from 
the surprising diversity of extrasolar planets, we should 
not assume that our irregular satellite complement is typ- 
ical, or that extrasolar analogues should follow all the 
same rules. 

Thankfully, some of the most important irregular 
satellite properties are dynamical and would have been 
discovered even if the Solar System had no irregulars. 
The main dynamical curiosity is their inclinations, which 
are all within about 60° of the ecliptic (but include retro- 
grade orbits). This evacuation of near-polar orbits is due 
to Solar and planetary perturbations, which drive the ec- 
centricities of highly inclined orbits to such large values 
that they either encounter regular satellites or leave the 
Hill sphere (Carruba et al. 2002; Nesvorny et al. 2003). 

Another constraint comes from the stability of cir- 
cumplanetary orbits. Although Shen & Tremaine (2008) 
find that satellites out to a few Hill radii could sur- 
vive the age of the Solar System around Uranus and 
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Neptune, all currently known irregulars have orbits 
with semi-major axes a p i less than half the Hill radius 
Rh ( e.g. Sheppard & Jewitt 2003; Holman et al. 2004; 
Sheppard et al. 2005) 

Rk = a p i(M pl /3M*) 1/3 , (1) 

where M p i is the planet mass and M* the stellar mass. 
On the sky, the giant planets' Hill radii span several 
degrees. The exact stability limit has a small inclina- 
tion dependence in that retrograde orbits are stable at 
larger distances than the widest stable prograde orbits 
(Nesvorny et al. 2003). 

In this section we outline a model for the colli- 
sional evolution of irregular satellite swarms. Because we 
want to make predictions of the only possible extrasolar 
observable — dust — we keep our model simple. There are 
many uncertainties in extrapolating a swarm of irregular 
satellites to a cloud of dust, such as the strength of satel- 
lites, the size distribution slope, and the minimum grain 
size, which at this stage make the development of a more 
complex collisional model largely unnecessary. 

The next four subsections contain many equations 
that describe properties of a circumplanetary swarm. 
Readers looking for actual numbers may like to refer 
ahead to Table 1, which shows estimates of some proper- 
ties for the Solar System giant planets. 

2.1 Collisional mass loss 

In a steady-state collisional cascade the mass within a 
given size range decreases as these objects are destroyed 
in collisions, but is replaced at the same rate by fragments 
created by destruction of larger objects. Mass is lost at 
the bottom end of cascade, usually by radiation forces 
that remove grains smaller than some minimum size. The 
evolution of the size distribution is therefore dictated by 
the collisional decay of the largest objects. 

The result of such a collisional cascade would, in an 
ideal situation (cascade infinite in extent, strength inde- 
pendent of size) , have a steady state size distribution with 
a well defined slope of n{D) = K D 2 ~ 3q where n(D)dD 
is the number of satellites between D and D + dD and 
q = 11/6 (Dohnanyi 1969). 

In fact the true distribution of circumstellar colli- 
sional cascades like the asteroid belt is expected to have 
different slopes in different size ranges due to the way 
strength depends on size (O'Brien & Greenberg 2003). 
Strength is typically described by the parameter Qd, 
which is the kinetic collision energy per target mass 
needed to shatter and disperse the target, such that the 
largest remnant is half the mass of the original target. 
Such a collision is commonly termed "catastrophic." 

Small objects are held together by their own mate- 
rial strength, and grow weaker with increasing size due 
to the increased likelihood of the presence of a significant 
flaw. To quote Benz & Asphaug (1994), "Subdivide this 
same rock into 100 equal pieces and 99 of them are now 
stronger than the original, owing to the simple fact that 
they do not contain the one weakest flaw." Above some 
transition size (D t ~ 0.1km) bodies gain strength from 
self-gravity. The energy needed for catastrophic disrup- 



tion now increases with size. Though objects may be shat- 
tered, extra energy is required to ensure the fragments 
have sufficient escape speeds and are no longer bound. 
Gravity also limits fracture propagation within the ma- 
terial, thus adding strength (e.g. Benz & Asphaug 1999). 
This behaviour is usually modelled using complex numer- 
ical codes, and parameterised by a power law for each of 
the strength and gravity regimes. In fact varies by 
about a factor 10 over the range of impact parameters 
and there are similar differences between strength laws 
derived by different studies (e.g. Benz & Asphaug 1999; 
Stewart & Leinhardt 2009). Thus Q^, and the resulting 
size distribution are the most uncertain inputs for our 
model. 

We set planetesimal strength with the 
Benz & Asphaug (1999) law for ice at 3 km s _1 . 
For objects larger than D t — 0.1 (in km), 

Ql = 0.1 pD 126 /f Q (2) 

in J kg -1 where the mass density p is in kg m~ 3 . The 
strength dependence of small objects (ocD~ ' 39 ) is only 
used in setting the size distribution of objects smaller 
than D t . Following Bottke et al. (2010), we allow objects 
to be weaker than the Benz & Asphaug (1999) law by in- 
cluding the factor /q (see also Levison et al. 2009). The 
strength law is similar to the Stewart & Leinhardt (2009) 
strength law when /q = 8. Bottke et al. (2010) repro- 
duce the Solar System's irregular satellite populations 
best when /q > 3 so we set /q = 5. 

The size distribution of objects with such strength 
properties is expected to have a slope with g s = 1.9 
at D < Dt (in the strength-dominated regime) and 
q g = 1.7 for D > D t (in the gravity-dominated 
regime, O'Brien & Greenberg 2003). Although sev- 
eral wiggles are also expected in the distribution 
(Campo Bagatin et al. 1994; Durda et al. 1998), here the 
size distribution is assumed to be continuous with the ap- 
propriate slopes (q s or g g ) between the smallest objects 
of size D m in (in pm) and the largest objects participating 
in the collisional cascade of size D c (in km). As we show 
later, D c may be smaller than the largest object, which 
has size -D max . 

With this two phase size distribution, the conversion 
between the size distribution's surface area (cr tot in AU 2 ) 
and mass (M to t in Mg) is 

Mtot = 0.0025 patot Dt~ 3q - D?*" 3 * D^ 5 (3) 

D 3G/ g 

where we assume D m i n < A < D c . For q s = 1.9, q s = 
1.7, and D t = 0.1 (in km), this relation simplifies to 

M tot = 3.9 x 10" 6 p atot D° c 9 D°J m . (4) 

This mass only includes objects between Z? m i n and D c 
(not D m ax) because we use M to t below to calculate colli- 
sion rates. 

The collisional lifetime of satellites of size D c can be 
calculated from the total mass using the particle-in-a-box 
approach. Here we follow Wyatt et al. (2010), who also 
take into account that objects of size D c can be destroyed 
in catastrophic collisions by those down to a size X c D c 
(assumed to be 3> 0.1km) where X c = (2Qo/v 2 el ) 1//3 
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and the collision velocity w re i is in m s . The rate of 
catastrophic collisions is 

6-3g g u re i CiX c ° 2 Mtot 



It 



8.4 x 10" 



3<& 



(5) 



in years 1 , where V is the volume occupied by the satel- 
lites in AU 3 . The inverse of the collision rate is called the 
collision time t c . The assumption of Dt <C -Xc D c means 
that only g g is needed for the collision rate. To obtain 
equation (5), it was necessary to integrate over the size 
distribution from X c D c to D c , which yields a function 
called G{q,X c ) by Wyatt et al. (2007). In Wyatt et al. 
(2010) the approximation G(11/6,X C ) « CiAf 2 with 
Ci = 0.2 and C2 = —2.5 is used (the limit of small X c ). 
Substituting these values and q s = 11/6 yields Equa- 
tions (9) & (10) of Wyatt et al. (2010). Here, we find 
a numerical approximation for the function G(q,X c ), 
which is within 10% for the more physically plausible 
range 0.01 < X c < 0.75 and 1.7 < q < 2, yielding 
d = 2.62(9-1.66) andC* 2 = 2.70(0.98-g). For q g = 1.7, 
Ci = 0.1 and C 2 = -1.9. 

Irregular satellites are assumed to orbit the planet at 
semi-major axes relative to the Hill radius in the range 
77 ± dri/2 (we use drj — 77/2) . The volume the satellites 
occupy is V = 47r?)s d^srjl x 0.866, where the extra factor 
of 0.866 accounts for the lack of near polar orbits. 

The mean relative velocity of collisions is expected 
to be some fraction / Vrol of the Keplerian velocity at 77, 
which (in m s _1 ) is 



(6) 



w k = 516M p 1 1 /3 Mi /6 /(r ? a p i) 1/2 . 

That fraction will depend on the eccentricities and in- 
clinations of the satellite swarm. A simple estimate of 
the mean collision velocity comes from assuming circu- 
lar and isotropic orbits, yielding / VraJ = 4/-7T and typ- 
ical impact velocities of ~0.5-3km/s for Solar System 
giant planets. These velocities are high enough that the 
impactor/target mass ratio for catastrophic collisions is 
small, so the energy lost in a collision is also small. Thus 
very little kinetic energy is lost in a typical collision and 
collisional damping is unimportant. Using the Monte- 
Carlo eccentric ring model of Wyatt et al. (2010), we find 
that the mean collision velocity is similar for a realistic 
orbital distribution, and somewhat lower when eccentric- 
ities are introduced. The values vary between about 0.9- 
1.3, so we adopt / Vrol = 4/n throughout. In fact orbits 
have a range of inclinations and eccentricities, and each 
collision has a different probability, which is itself a func- 
tion of the relative velocity (Bottke et al. 1994). 

Substituting our approximation for C\ and C2 for 
q g = 1.7 yields the rate of catastrophic collisions 



Rc 



1.3 x 10 7 



(7) 



in years -1 . As one expects, the rate is independent of a p i 
for the same physical swarm (i.e doubling a p i halves 77). 
The factors that largely set the collision rate are 77, D c , 
dpi , and / Vrol . The strongest contributions are from 77 and 
dpi, which set the cloud volume and space density. The 
rate also depends strongly on the mean collision velocity 
because this speed sets both the rate at which an ob- 



ject sweeps through space and the number of impactors 
that result in a catastrophic disruption. Greater collision 
velocities mean smaller impactors can destroy a given ob- 
ject, and smaller impactors are more numerous. Because 
Qd also depends on D c (eq. 2), the collision rate also de- 
pends strongly on D c . For the same total mass, larger D c 
means fewer large objects, which are also stronger. 

Curiously, the least important parameter is the mass 
of the planet itself. This result arises because in equa- 
tion (7), ?j rc i is slightly less than linearly dependent on 
M p i, which nearly cancels with the linear dependence of 
volume on M p \. In simple terms, for fixed 77 and M to t 
the increase in volume with planet mass works out to be 
slightly stronger than the increase in velocity, so the col- 
lision rate decreases slowly as the planet mass increases. 

We validate our model with the eccentric ring model 
(Wyatt et al. 2010). Our collision rate is within a factor 
of three for a range of eccentricities with the best agree- 
ment for high values, sufficient for our purposes here con- 
sidering that much larger uncertainty lies with assumed 
material properties and the resulting size distribution. 



2.2 Radiation forces on dust 

To derive the surface area in small grains (dominated by 
small objects) from the total mass (dominated by large 
objects) we need to know the size of the smallest grains 
that can survive in circumplanetary orbits. Two notable 
detections are micron-sized grains found orbiting at large 
(50-350 .Rjup) distances from Jupiter (Krivov et al. 2002; 
Kriiger et al. 2010), and the large tenuous ring found 
orbiting Saturn (Verbiscer et al. 2009). Both studies at- 
tribute material released by impacts from interplanetary 
grains as a likely source, though Verbiscer et al. note that 
debris from irregular satellite collisions impacting Phoebe 
could also be the cause. 

As with grains orbiting a star, the effect of radiation 
forces on dust characterised by /3 = irradiation /-Pgravity 
(both due to the star) plays the most important role in 
setting the minimum size of grains that survive on cir- 
cumplanetary orbits (Burns et al. 1979). Other effects re- 
lated to interaction with planetary magnetospheres (e.g. 
Horanyi 1996) play some part but are less important for 
the >/-im grain sizes and wide orbits considered here. 



2.2.1 Radiation pressure 

Radiation pressure is the radial component of the force, 
which in contrast to circumstellar orbits causes the or- 
bits of dust grains to evolve. While semi-major axes re- 
main constant, eccentricities oscillate with a period equal 
to the planet's orbital period, with a maximum that de- 
pends on /3 and the grain orbit. The maximum j3 of grains 
that survive in orbit around the planet with e < 1 have 
(Burns et al. 1979) 



/3c = v/3v Q = 5.8 x 10" 



3 m p Y 3 



/(M 



1/3^1/2) 



(8) 



where v is the velocity of a grain as it orbits the planet, 
and Vq is the velocity of the planet as it orbits the star. 
For typical planets and irregular satellite orbits /3 C is 
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much smaller than the blowout limit of 0.5 for stellocen- 
tric orbits. 

Because ft for normal grains peaks where the star 
radiates most of its radiation, it might be possible for /3 C 
to allow both large and very small grains to survive, with 
only grains in the peak being excluded. However, for sub- 
micron "astronomical silicate" grains, those on the small 
side of the peak, the smallest grains have j3 ~ 0.11 (e.g. 
Gustafson 1994) and is higher for more massive stars, 
so grains smaller than the wavelength of typical stellar 
radiation will usually be ejected. 

To convert p c into a minimum size, we use D m in — 
(1150/p/3 c )WM* in fj,m (e.g. Wyatt 2008a). The mini- 
mum size is therefore 



D min = 2 x 10" 



n 



1/2 



p Mil 3 Ml' 3 



(9) 



in fim.. For typical parameters, Z) m m is at least /im-size. 
For Jupiter and Neptune, Equation (9) yields 12 and 
23/im, an order of magnitude larger than the minimum 
(blowout) size for the same grains on circumsolar orbits. 



2.2.2 Poynting-Robertson drag 

An alternative to grain removal by radiation pressure is 
orbital decay due to Poynting-Robertson (PR) drag. The 
decay timescale is similar to the heliocentric case 



t PR = 530ap 1 /(/3M* 



(10) 



in years (Burns et al. 1979). However, as noted above /3 
for the smallest grains is typically much smaller than 0.5 
and tpn correspondingly longer. 

For grains to spiral into the planet by PR drag 
they must avoid colliding with other grains first, which 
breaks them into smaller particles that are instead re- 
moved by radiation pressure. The competition between 
PR drag and collisions can be characterised by xpr = 
f pr /trad- Wyatt et al. (1999) showed that for the small- 
est grains this collision rate is roughly t peT r^ ust /(4: cr tot ) 
for a flat disk with radial extent rd us t ± J"d us t/4 (where 
t pcr is orbital period). Adapting this expression to an 
isotropic case results in a small change due to the greater 
cloud volume and faster collision velocities: t C oi,dust = 
t peI rf vlSt /(4 f Vxel (7tot), or using our parameters 



= 10" 



(ya pl ) 7/2 M 2 J 3 
M 7 J 6 a tot 



(11) 



in years. Expressed in terms of our basic parameters, the 
ratio is 



7/6 



Xph 



4 x 10 



4 P Dmin Ctot M, 

a 3 J 2 V 7/2 M 2 J 3 L+ 



(12) 



When this ratio is larger than unity, grains suffer col- 
lisions before their orbits have time to decay due to PR 
drag, and are subsequently removed by radiation pressure 
(termed "collision dominated"). When this ratio drops 
below unity, grains spiral in towards the planet before 
they collide and may encounter any existing regular satel- 
lites as they do so ("PR dominated"). Substituting Equa- 
tion (9) for Dmin yields 



Xpr. = 8 x 10 



, 9 a tot Ml /2 



a 3 J 2 T] 3 M p i 



(13) 



With only the Solar System example to go by (see §3.1.1), 
it is impossible to tell whether the "typical" satellite dust 
cloud will be collision or PR dominated. However, due to 
the amount of dust needed for an extrasolar irregular 
satellite swarm to be detectable (crtot > 10 -4 AU 2 ), any 
observed extrasolar swarms will likely be collision domi- 
nated. 



2. 2. 3 Summary 

In order to estimate the minimum grain size we have nec- 
essarily made a number of simplifications. We have used 
expressions for radiation forces assuming low inclination 
orbits akin to Saturn's ring. In fact /3 C varies with the ori- 
entation of the orbit relative to the Solar direction and 
is smallest for those with pericenter initially aligned with 
the Solar direction. Equation (8) assumes a coplanar or- 
bit with pericenter perpendicular to the Solar direction 
and underestimates /3 C for certain loss by a factor 2-3 
(Burns et al. 1979). We therefore overestimate the mini- 
mum grain size, and underestimate a to t- 

Our prescription for Z) m m does not preclude detec- 
tion of grains smaller than D m i n . Grains below the min- 
imum size may complete many orbits before their ec- 
centricity exceeds unity. Also, the minimum grain size 
decreases with r\ so small bound eccentric grains with 
smaller planetocentric semi-major axes can be found any- 
where within 2rj. On a detailed level, _D m i n takes on a 
range of values and for the smallest grains with /? ~ f3 c , 
the orbit averaging used in deriving Equation (8) breaks 
down. Using numerical simulations, Krivov et al. (2002) 
find that D m i n for Jupiter is ~1 fim (compared to our 
value of 12.4 /im, see Table 1). Considering that our value 
is overestimated by the factor 2-3 noted above, the min- 
imum size is probably 1 to a few /im, with differences in 
the assumed grain properties contributing some uncer- 
tainty. 

While our simplifications are reasonable, they gloss 
over important aspects of grain dynamics. For example, 
Krivov et al. (2002) show circumplanetary dust clouds to 
have both size-dependent and pro/retrograde orbit sen- 
sitive structure. Because non-gravitational forces cause 
small objects to deviate from the orbits of their parent 
bodies, such forces lead to effects beyond the scope of our 
model. Therefore, future work will need to consider how 
grain orbital evolution affects both the physical appear- 
ance of the cloud, and the underlying size distribution. 
For example, in a typical collision dominated circumstel- 
lar disk the minimum (blowout) grain size is a single 
number, independent of radial distance. This fairly sharp 
truncation creates a wave in the size distribution due to 
the alternating lack and then over-abundance of projec- 
tiles that destroy larger targets (Campo Bagatin et al. 
1994). However, in a circumplanetary cloud the mini- 
mum grain size varies with circumplanetary distance (Eq. 
9), which could lead to qualitatively different and spa- 
tially varying size distributions. In §3.1.2 we suggest that 
Jupiter's swarm may be PR dominated, which can lead 
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to further differences at the small end of the size distribu- 
tion. In addition, other non-gravitational forces known to 
be important for circumstellar dynamics (e.g. Yarkovsky 
effect, Burns et al. 1979) may be important for the evolu- 
tion of circumplanetary orbits and consequent collisional 
evolution. 



2.3 Time evolution 

We now turn to the time evolution of disk properties. 
Assuming that the size distribution remains fixed and 
that mass is lost by catastrophic disruption of the largest 
objects, the disk mass remaining as a function of time is 
found by solving (e.g. Wyatt et al. 2007) 



-^Opt 



= 50 



M 



0.33 r0.55 

J v i cl 



dMtot 



dt 

which yields 
M to t(i) 



— — Aftot Rc 



Mtot(0) 



l + i?cc(0)t 



(14) 



(15) 



where we use t or to indicate variables at a particular 
time where necessary. The initial cloud mass is Af to t(0) 
and Mtot(t) = Mtot(0)/i? cc (0) t when R cc (0)t > 1. 
Equation (4) shows that mass and surface area are lin- 
early proportional for fixed D c and -D m in so the surface 
area of small grains decays in the same way. Because the 
collision rate depends on the remaining mass, after a few 
collision times the remaining mass is independent of the 
initial mass (Wyatt et al. 2007). 



2.3.1 Where are the most massive swarms? 

There is an interesting interplay between the initial cloud 
mass and planet semi-major axis. For arbitrarily large ini- 
tial masses, swarms at any distance and time are on the 
1/t part of their evolution where the remaining mass is 
independent of the initial mass (called a "collision lim- 
ited" disk by Heng & Tremainc 2010). Because these sys- 
tems have a maximum remaining mass that depends on 
the collision time (l/7? cc ) without the mass term, the re- 
maining mass increases strongly with planet semi-major 
axis. In fact, t c oc a^ 13 (i.e. more strongly than a^), so 
collision limited swarms are brighter in scattered light for 
larger planet semi-major axes (for fixed planet mass, age, 
etc.). 

Of course, the initial swarm mass cannot be arbi- 
trarily large, so swarms around sufficiently distant plan- 
ets will take some time to start colliding. Thus, for fixed 
initial mass, swarms around close planets will rapidly de- 
cay due to the short collision times, while more distant 
planets all have the same mass in satellites because the 
largest irregulars have not yet, or only just begun to col- 
lide. Of these more distant planets, the one whose swarm 
has just started to suffer collisions is the one that receives 
the most stellar insolation and is brightest in scattered 
light. 

The semi-major axis of the planet hosting this swarm 
can easily be worked out from the collision rate in equa- 
tion (7), because a swarm that has just started to collide 
has t = 1/R CC , therefore 



n 



tMtot(O) 
pD c 



(16) 



in AU. The "opt" subscript indicates that this planet 
lies at the optimum distance to be detected at "opt"-ical 
wavelengths. This prescription for the brightest swarm 
is complicated by the fact that faint objects are harder 
to detect close to host stars, an issue we return to in 
§3.2. The same concept of an optimal distance applies to 
thermal emission, but is somewhat different because the 
cloud temperature changes with planet semi-major axis. 

2.3.2 Stranding the largest objects 

If collisional evolution proceeds for long enough, the mass 
in the swarm will drop to the point where it is com- 
parable with the mass contained in a single largest ob- 
ject. Around this time the evolution of the largest ob- 
jects changes from being reasonably well described by our 
particle-in-a-box formalism, to a regime where individ- 
ual collisions and cratering are important (Bottke et al. 
2010). In this regime, the largest objects are less likely to 
be destroyed due to their small number and the decreased 
number of potential destructors. While our model cannot 
take cratering or stochastic collisions into account, we can 
approximate the evolution by assuming that when the 
number of largest objects drops too low, they lose their 
connection with the rest of the size distribution and be- 
come "stranded." The size of the largest non-stranded ob- 
ject is of size D c by definition, which decreases over time 
and leaves a relatively flat size distribution of stranded 
objects between sizes D c and D max , whose evolution is 
halted due to a lack of would-be destructors. 

We implement this simple approximation by assum- 
ing an object of size D is stranded when the number of 
objects between sizes D/2 and D 



i(D/2 -> D,t) = 



K (2 3 ^- 3 _ i) (io 3 Df~ 3q « 
3<? g - 3 



(17) 



drops to some number A r str . This number sets the normal- 
isation for the size distribution of stranded objects. This 
assumption sets the size distribution slope of stranded 
objects g str = 1 because n(D/2 — > D,t) is independent 
of D for D c > D > D 

max ■ 

The first object is stranded at f n icft, which can be 
calculated from the initial number of objects in this size 
range 



n(D c /2 -+D c ,t = 0) 

Rcc(0)N str 



(18) 



The difference between t co i and t n i c ft is simply a measure 
of the initial number n of large objects, which take n/N s t T 
collision times to become stranded. Like the remaining 
mass once collisions occur, this time is independent of 
the initial cloud mass. 

At this point, D c and Dmax become distinct sizes, 
with Dmax remaining fixed and D c decreasing with time 
as smaller objects are stranded. The remaining planetes- 
imal population decays at the collision rate for D c size 
and strength objects and the mass remaining in the size 
distribution below D c . Therefore, by substituting D c for 
Mtot in equation (14), the evolution after t n ieft obeys 
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Figure 1. Example evolution of the size distribution (left) and total surface area cr tot (right) for lOGyr. In the left panel the lines 
become lighter for later times and are logarithmically spaced in time. The dashed line shows stranded objects. In the right panel 
the evolution with stranding is shown as a solid line and the "normal" 1/t evolution with no stranding as a dashed line. 



dt 



-Dt R c 



(19) 



where we have used Mtot oc Dj? (for fixed n(D c /2 — > D c )). 
The collision rate is for size D c objects, on which both 
A/tot and Qd depend, and varies as R cc oc D\ 2 (eqs. 2 & 
7). Integrating equation (19) yields 



(1 + 0.4(f - tnleft)/£„left) a 



(20) 



where a = 1/1.2. This evolution is illustrated in Figure 
1, which shows the time evolution of the top section of 
the size distribution and crtot. The size distribution ini- 
tially decays straight down (i.e. K decreases) with D c 
fixed. When there are only AW largest objects left, _D C 
begins to decay as dictated by Equation (20). The size 
distribution then moves to the left (smaller D c ), with 
both crtot and M to t continuing to decay. The right panel 
shows that crtot drops more quickly after the first objects 
are stranded, tending to cr to t oc t -1 ' 75 (using eqs. 4 and 
20). Such evolution is potentially interesting, as a dust 
cloud is accelerated towards being PR-dominated after 
stranding due to the stronger decrease in small grains 
with time. Entering the PR-dominated regime, the size 
distribution is effectively truncated at larger sizes, lead- 
ing to an even faster decay of cr to t than shown in Figure 
1 (Dominik & Decin 2003). 

This evolution is necessarily very simple because ob- 
jects are only destroyed by catastrophic collisions in our 
model. However, at these late stages the mass released 
into the cascade by cratering may be as or more impor- 
tant (Bottke et al. 2010). Therefore, though our model 
is physically plausible, the actual evolution will depend 
on details such as the relative importance of catastrophic 
disruptions vs. cratering or differences between pro and 
retrograde populations. We treat N s t T and a as free pa- 
rameters when comparing our model with the Solar Sys- 
tem irregulars in §3.1. 



2.4 Observable quantities 

Given the surface area in dust, we derive the flux den- 
sity F due to the cloud and planet from both scattered 
light and thermal emission. Because distance and surface 
area can be in different units (e.g. m, pc, AU), quantities 
in these equations have dimensions. We take Solar Sys- 
tem planetary effective temperatures and radii from Cox 
(2000). The stellar flux at the planet is 



F+ = i* B4\,n)/(4,a K T?al 



(21) 



where <tk is the Stefan-Boltzmann constant. The scat- 
tered light seen from Earth is (e.g. Collier Cameron et al. 
2002) 



F scat = F+AgQ/iirdli) 



(22) 



where A is either the projected area of the planet or tj to t 
for dust. The geometric albedo Q is assumed to be 0.08 
for dust, similar to both the values for Kuiper belt ob- 
jects, Jovian trojans, and irregular satellites in the So- 
lar System (Stansberry et al. 2008; Mueller et al. 2008; 
Fernandez et al. 2009) and to that inferred for the Fo- 
malhaut dust ring (Kalas et al. 2005). For planets we 
use an albedo of 0.5. The phase function g is set to 
unity for the Solar System because planets are on ex- 
terior orbits. For extrasolar planets and swarms we set 
g — 0.32, the value for a Lambert sphere at maximum 
extension from the host star (e.g. Collier Cameron et al. 
2002). We use a blackbody estimate for the dust tem- 
perature Tdu S t = 278.3 Li /a*; K. The distance d p i is 
the distance to the system in the extrasolar case, or a p \ 
in the Solar System (i.e. roughly the distance for a Sun- 
Earth-Planet angle of 90° for outer planets) . For thermal 



F t h = B v {\, T) A/d 2 pl . 



(23) 



For dust we include a "greybody" decrease in emission of 
210/A beyond 210 to account for inefficient emission 
from small grains (e.g. Wyatt 2008b). Though the actual 
spectrum depends on the grain properties and size dis- 
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tribution, this addition provides more realistic (sub-)mm 
flux densities than a plain blackbody. 



3 APPLICATIONS 

In this section we apply our model to three different ir- 
regular satellite populations. We first compare our model 
evolution with the Solar System's complement of irregu- 
lars and make predictions of the current levels of dust. We 
suggest that the fate of dust changes over time, with im- 
plications for regular satellites such as Callisto and Iape- 
tus. We then apply our model to possible satellite clouds 
around extrasolar planets. Finally, we apply our model 
to Fomalhaut b. 



3.1 Solar System Irregular Satellites 

In this section we compare our model with the Solar 
System's irregular satellites. We model a swarm with 
initial mass M to t(0) = 0.01M moon . This mass ensures 
all swarms are collision limited and therefore does not 
influence the model. We set D c = 150 km for Jupiter 
and Uranus and D c — 250 km for Saturn and Neptune. 
We set the density similar to known irregular satellites 
p = 1500 kg m" 3 (e.g. Jewitt & Haghighipour 2007). 
Though we assume f Vxel — 4/n at each planet, this num- 
ber varies somewhat depending on the specific orbital 
properties of each swarm. We use the average n of known 
irregulars at each planet. 

Figure 2 shows a comparison of our model with the 
Solar System's irregular satellites. 1 We have combined 
the pro/retrograde satellites for this distribution. The 
simplicity of our model means that we cannot account 
for the different evolution of pro/retrograde populations, 
which appears to be important at Uranus (Bottke et al. 
2010). Because all satellites at each planet are at a sim- 
ilar distance from Earth, these distributions are near- 
complete to the smallest observed size and need no cor- 
rection (e.g. Sheppard et al. 2005). 

We vary the normalisation iV str and rate of decay 
after stranding power-law index a to obtain a by-eye fit. 
Variation of these parameters over a fairly small range al- 
lows an excellent match for Jupiter, Saturn, and Uranus. 
Given the differences in the irregular populations at each 
planet and the simplicity of the prescription for strand- 
ing, some variation is expected. The first few columns 
of Table 1 show the rj, t co i, and t n \ e ft for each planet. 
The slower collisional evolution for more distant planets 
means that t n i e ft is longer and that the size distribution at 
the current epoch turns up (i.e is stranded) at larger sizes. 
The minimum known satellite size increases for more dis- 
tant planets, making the model comparison less certain 
for Uranus and Neptune. We do not know if their size 
distributions are similar to Jupiter and relatively flat to 
~8 km, or instead turn up at larger sizes like our model. If 
evolution at Uranus is similar to Jupiter and Saturn, our 

1 Taken from http:/ /www. dtm.ciw.edu/users/shcppard/satellites/ 
in March 2010 



model predicts that the size distribution should turn up 
at a few tens of km due to the slower collisional evolution. 

There are few Neptune irregulars to compare 
with, perhaps because many were depleted by Nereid 
(Nesvorny et al. 2003) and/or Triton (Goldreich et al. 
1989; Cuk & Gladman 2005). Nereid's low orbital incli- 
nation (7° relative to the ecliptic) and relatively close 
prograde orbit (n = 0.05), as well as a colour and albedo 
similar to the Uranian satellites Umbriel and Oberon 
(Buratti et al. 1997), mean that it may in fact have 
formed as a regular satellite of Neptune. For these reasons 
Bottke et al. (2010) did not model Neptunian irregulars. 

Given these complications, it is perhaps unsurpris- 
ing that our model for Neptune's irregulars needs very 
different N s tr and a to the other giant planets, and is 
still a poor match. This difference is due to the slow 
collisional evolution, which predicts that t n i e ft is a size- 
able fraction of the Solar System age. With stranding 
occurring at such a late time, the subsequent evolution 
must be very rapid to deplete the population to that cur- 
rently observed. With values for N stI and a more like 
the other three planets, our model predicts that Neptune 
would have two orders of magnitude more satellites (i.e. 
be at about the second to lowest curve) . This discrepancy 
suggests that either the initial conditions for Neptune's 
swarm were quite different to the other planets, or as 
already proposed, the irregulars were affected by Nereid 
and Triton. 

Based on these comparisons we conclude that our 
model provides a reasonable description of the collisional 
evolution of an irregular satellite swarm. 



3.1.1 Current irregular dust levels 

Given the ongoing collisional erosion of the irregular 
satellites, the presence of dust is inevitable. Here we 
make some estimates of the expected level of dust at 
each planet and relate them to a few relevant observa- 
tions. Table 1 shows estimates for each planet based on 
<7tot from the models of Figure 2 at t = 4.5 x 10 9 years. 
That is, they are extrapolated using our size distribution 
and independent of the collisional model. We calculate 
the flux density and surface brightness at 1 and 100 /im 
(distributed uniformly over a disk of radius 0.5-Rh for 
simplicity), which roughly correspond to peaks in scat- 
tered light and thermal emission respectively. 

The predicted surface brightness levels are much 
fainter than the background in the ecliptic, but expected 
to vary on a similar scale (i.e. degrees). For comparison, 
the Zodiacal background at 1.25 and 100 fim is about 0.4 
and 9MJy/sr respectively (Kelsall ct al. 1998). Detection 
of these clouds therefore requires accurate subtraction of 
this (and the galactic and cosmic) background. 

Of course dust that originates from irregular satel- 
lites has been detected already (see also §3.1.2 below). 
The largest of Saturn's rings is probably fed by mate- 
rial generated when inter/circumplanetary grains impact 
the irregular satellite Phoebe (Verbiscer et al. 2009). At 
24 /xm, this ring has a surface brightness of ~0.4MJy 
sr _1 , which is less than 1% of the (~70 MJy/sr) zodiacal 
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Figure 2. Comparison of our model (thick curves) with the Solar System's irregular satellites (dots and thin lines, pro/retrograde 
populations have been added together). Curves show ten logarithmically spaced times between 10 6 and 4.5 X 10 9 years (dark to 
light). The legend indicates the planet and values for 7V s t r and a. Model curves turn down at the largest sizes because they are 
derived from integrating n(D) (and there are zero objects larger than -D ma x). 
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0.4 


2.5 
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12 


1.5 


0.9 


0.05 


12.0 


37000 


0.0034 


960 


450000 


0.26 


Saturn 


0.43 


0.3 


17.0 


740 


16 


7.5 


1.5 


0.8 


5.7 


2400 


0.0035 


840 


60000 


0.52 


Uranus 


0.46 


0.2 


14.0 


1100 


21 


16.0 


1.1 


13.0 


0.8 


28 


0.0016 


210 


950 


0.45 


Neptune 


0.77 


0.2 


250.0 


3500 


23 


36.0 


0.6 


13.0 


0.3 


■1 


0.0005 


99 


350 


0.19 



Table 1. Irregular satellite model and dust properties at each planet. Estimates are based on the models in Figure 2. The last six 
columns show estimates of the total dust cloud and planet flux density at opposition, and the surface brightness if the cloud were 
evenly spread over a disk with radius 0.5-Rhj a t 1 and 100 fim. The peak surface brightness is ~6 times higher than shown here if 
grains are distributed as in Figure 3. 



background. For comparison, our Saturn swarm model 
has ~0.2 MJy sr _1 when spread over a 0.5 Rh radius disk. 

To take a more detailed look at Saturn, Figure 3 
shows an example of what such a cloud might look like 
from Earth at 24 (im. In creating this image we have as- 
sumed that the dust follows the same orbits as the par- 
ent bodies, which have semi-major axes distributed be- 
tween 0.04-0.16 AU, and e = 0.1-0.6 based on Figure 
7 of Nesvorny et al. (2007). The image would therefore 



look the same at any wavelength, though the total flux 
changes. The total flux in this 24 /im image is 320 Jy, 
which yields a peak surface brightness of around 1.25 
MJy/sr. Therefore, the surface brightness predictions in 
Table 1 are a factor ~6 too low if grains follow their 
parent body orbits. The peak level is about three times 
brighter than Saturn's Phoebe ring. 

Though the cloud is spread over a much larger region 
of sky, for it to have evaded detection thus far, particu- 
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Figure 3. Simple model of Saturn's irregular satellite dust 
cloud at 24 (im. Grains follow parent body orbits. The total 
flux is 320 Jy and the scale is a linear stretch between and 
1.25 MJy/sr. Flux from Saturn itself is not included. 



larly in the Phoebe ring observations, it appears that 
our prediction is at least a few times too high. However, 
Verbiscer et al. (2009) note an apparent trend of increas- 
ing surface brightness towards the Phoebe ring, which 
could be due to dust from irregular satellites. The nar- 
row vertical range of their Figure 3 does not constrain 
our prediction, because the Phoebe ring could be sitting 
on top of a larger background. 

Several uncertainties with our model may explain 
this apparently high prediction for the level of dust at 
Saturn, the most likely being that small departures from 
our assumed size distribution can lead to large differ- 
ences in the predicted surface area in dust. In addition, we 
have assumed that grains follow the orbits of their parent 
bodies, but Krivov et al. (2002) show that the detailed 
cloud structure is more complex. Finally, Table 1 shows 
that grains around Saturn may be in the PR-dominated 
regime, which can alter the cloud structure and reduce 
the dust level. 

Despite our apparent over-prediction of the cloud 
surface brightness, detection of such a cloud is difficult, 
not least because it may be unexpected. The large extent 
means achieving sufficient coverage (ideally the entire Hill 
sphere) is expensive at IR wavelengths (i.e. using space 
telescopes) . Poorly characterised dust bands will also hin- 
der background subtraction. At longer wavelengths, such 
as those covered by Herschel SPIRE, achieving decent 
coverage is easier, but the background comprises a com- 
bination of Zodiacal and galactic light. Care would be 
needed to ensure good background subtraction that min- 
imises elongation effects. 

Now looking at Jupiter, from Galileo dust detec- 
tion data Krivov et al. (2002) derive a constant dust 
number density of ~10km~ 3 between 50 and 350-R.j up . 
Kriiger et al. (2010) report detection of micron-size parti- 
cles (and note the lack of detections of smaller particles) 
by the Galileo dust detector at a distance of approxi- 



mately 350i?jup. Such large planetocentric distances are 
prime irregular satellite territory, being about half the 
Hill radius. These grains have previously been explained 
as ejecta from impacts of interplanetary dust on irregular 
satellites, which produce about the same level of dust as 
detected by Galileo (Krivov et al. 2002). The predicted 
space density of dust for Jupiter shown in Table 1 (again 
within 0.5 -Rh) is at a lower level to the Galileo detec- 
tions and may not have contributed to the measurements. 
However, the space density is an extremely strong func- 
tion of D m i n . For example, if the minimum grain size 
in our model were 1 /jm the predicted space density at 
Jupiter would be ~100km~ 3 . Therefore, Galileo obser- 
vations constrain either the grain size in our model to be 
larger than a few microns, or the dust level to be lower 
than predicted if the minimum grain size is 1 /im. 



3.1.2 Fate of Dust 

The best way to probe the small end of the irregular 
satellite size distribution is to detect the dust cloud di- 
rectly. However, as noted above such an observation is 
difficult and there are complementary ways to detect ir- 
regular dust. One of the most interesting signatures ex- 
ists on the surfaces of some regular satellites. The orbits 
of the smallest grains are strongly affected by radiation 
forces and may end up on a regular satellite. However, be- 
cause this deposition is only one of several possible fates 
a grain may meet, an understanding of which is more 
likely (and when) is needed to make a strong connection 
between irregular dust and regular satellite surfaces. Our 
model is too simple to model irregular satellite evolution 
in the Solar System at a detailed level. However, we offer 
some order of magnitude arguments concerning the fate 
of dust, which highlight questions that should be asked 
by more detailed studies. 

Like a circumstellar disk, a young circumplanetary 
swarm will most likely be collision dominated. At this 
stage grains are lost to the planet and interplanetary 
space. When the mass has been sufficiently depleted it 
becomes PR dominated and grains spiral in toward the 
planet. Grains destined to impact the planet, or pass 
nearby, may meet a third fate and be swept up by a reg- 
ular satellites. Though only a small fraction of mass may 
be lost this way, it is important due to the visible effect 
of leading/trailing asymmetries on some tidally locked 
Solar System regular satellites, such as Callisto and Ia- 
petus. These asymmetries are thought to arise from the 
higher accretion rate of retrograde dust by the leading 
hemisphere (e.g. Burns et al. 1996). 

Though most grains in collision dominated swarms 
either leave the Hill sphere or impact the planet due to 
radiation pressure, nearly all must pass through the reg- 
ular satellite domain to do so. Grains destined to hit 
the planet or leave the Hill sphere do not reach their 
maximum eccentricity immediately, but instead make a 
number of pericenter passages first. The smallest grains 
with r\ ~ 0.1-0.5 only complete a few, to a few tens 
of orbits before their eccentricities exceed unity (i.e. 



n 



3/2 



per .planet , 



and are therefore unlikely to 



encounter regular satellites before removal. 
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Grains in the small size range where e grows large 
enough to pass through the regular satellite region but 
not high enough to hit the planet or leave the Hill sphere 
have the best chance of colliding with regular satellites. 
However, this collision time must be shorter than the time 
for grains to collide with themselves. Taking Saturn as an 
example, at 77 = 0.3 with er tot = 7.5 x 1(T 9 AU 2 (Table 
1) the current time for collisions between the smallest 
grains is predicted to be ~10 6 years (eq. 11). In the past 
the level of dust was much higher, and the grain-grain 
collision time correspondingly shorter. 

Continuing with the Saturn example, grains with 
0.83 > e > 0.997 traverse the region between Iapetus' 
orbit and Saturn's surface near pericenter. This range 
corresponds to grains of sizes 16 > D > 25 /im. The or- 
bital period i per is about 3 years for these grains. Grains 
on coplanar orbits with Iapetus in this size range spend 
no more than 40% of their lifetimes in Iapetus-crossing 
orbits. Iapetus occupies about 1/20000 th of its orbital 
torus, so coplanar grains have a collision time of ~10 4 
orbits, or ~10 4 years. This time is shorter than the grain- 
grain collision time, so coplanar grains in this size range 
appear likely to impact Iapetus. 

However, grains on inclined orbits (i.e. the major- 
ity) only have a chance of hitting Iapetus if they cross 
Iapetus' plane at precisely the right radial distance. Ia- 
petus is ~1500 km across, compared to a semi-major axis 
of 3.5 x 10 6 km, so the chances of a plane-crossing par- 
ticle encountering Iapetus orbit are roughly 1/2000. The 
chance of an inclined grain impacting Iapetus at each 
plane crossing is ~10~ 7 per orbit, or an impact time of 
~10 7 years. The Iapetus collision time is therefore longer 
than the grain-grain collision time, so grains are more 
likely fragmented to smaller sizes first. Taking these num- 
bers as representative, it appears that in collision domi- 
nated swarms only a small fraction of grains, those within 
a small size range and on coplanar orbits, will impact reg- 
ular satellites. 

In contrast to the highly variable eccentricities of ra- 
diation pressure affected grains, PR drag causes grain or- 
bits to collapse slowly. The chance of impacting a regular 
satellite is larger than for the more distant and eccen- 
tric radiation pressure induced orbits because grains or- 
bit the planet more often as their semi-major axes shrink. 
For example, in the case of grains released from Phoebe, 
collision with Iapetus is nearly guaranteed, with escapees 
destined to hit Hyperion or Titan (e.g. Burns et al. 1996). 

Table 1 shows that dust from irregular satellites at 
Jupiter and Saturn is near or in the regime where PR 
drag dominates, but dust at Uranus and Neptune is not. 
If grains in PR dominated swarms are more likely to 
impact regular satellites, this finding is consistent with 
the brightness asymmetries seen on Callisto and Iapetus, 
which are less marked (but still present) on the Uranian 
satellites. 

While we suggest that the Jupiter and Saturn 
swarms are PR dominated now, they were collision dom- 
inated in the past. The possibility that collision domi- 
nated swarms do not coat regular satellites as efficiently 
as PR dominated ones therefore has implications for the 
interpretation of the brightness asymmetries. If the mass 
deposited on regular satellites is some fraction of the total 



mass lost, this fraction will increase as the swarm changes 
from collision to PR dominated. This change will affect 
the accretion history and more mass may be accreted at 
late stages if the difference in accretion efficiency is large 
enough. 

Such interpretations may be too simple. For exam- 
ple, Iapetus' asymmetry is likely due to accretion of 
grains from Saturn's Phoebe ring, so is due to grains re- 
leased from an individual irregular satellite rather than 
the cloud in general. In summary, further study of the 
fate of irregular satellite debris and potential observables 
requires consideration of the competition between colli- 
sions, PR drag, and radiation pressure. 



3.2 Irregular Satellite Clouds around 
Exoplanets 

While the dust produced by Solar System irregular satel- 
lites is at very low levels, these clouds were many times 
brighter at earlier times. In this section we explore the 
prospects for discovery of young extrasolar circumplan- 
etary swarms. Like circumstellar debris disks, dust can 
be discovered in scattered light at optical wavelengths, 
and thermal emission at IR wavelengths. In the case of 
thermal emission, one may need to distinguish a plane- 
tary atmosphere from dust at a similar temperature via 
spectral features. 

Because there is little a priori reason to choose any 
particular system configuration, we first model a par- 
ticular system and then show which configurations are 
detectable in a more general sense. Clearly, to be de- 
tectable with current or near-future technology these sys- 
tems must possess more dust than predicted for the Solar 
System's irregulars. Further, because the first objects dis- 
covered where sensitivity is an issue tend to be the most 
extreme (e.g. hot Jupiters), we do not restrict the ini- 
tial swarm masses. While the initial conditions used by 
Bottke et al. (2010) based on the Nesvorny et al. (2007) 
simulations were at most a small fraction of the Moon's 
mass, extrasolar systems may have much more massive 
planetesimal belts from which to capture satellites, and 
mechanisms disfavoured for the Solar System may also 
operate. We discuss capture mechanisms further in §3.2.1. 

In this section we first use the example of a Jupiter- 
mass planet orbiting an A5 star at 10 pc. We choose 
this spectral type of star simply because the only plan- 
ets directly imaged to date orbit A-stars (Kalas et al. 
2008; Marois et al. 2008; Lagrange et al. 2009). However, 
more massive stars do not necessarily produce brighter 
clouds because their high luminosity and mass increases 
Dmin and changes the collision rate. For the evolution 
of the planet, we use the non-irradiated Jupiter-mass 
(Z = 0.02) model from Baraffe et al. (2008). We set 
Aftot(O) = lOMmoon, and D c = 250 km. Another choice 
is when the satellite swarm is captured (or otherwise 
formed) and evolution begins. For simplicity, we assume 
swarms start evolving at t = 0. 

Figure 4 shows the evolution of scattered light and 
thermal emission from both planet and dust cloud when 
the planet orbits at 50 AU. The Figure is drawn to high- 
light where the spectrum of each component moves over 
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Figure 4. Example satellite swarm evolution around a 1 Mj up 
cxtrasolar planet at 50 AU around an A5 star between 10 7 
(darkest curves) and 10 10 years (lightest curves, logarithmi- 
cally spaced). A small part of the stellar spectrum is visible. 
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Figure 5. Adopted instrument contrast ratios. 



time (which is downwards in this plot). In scattered light, 
the dust cloud is initially much brighter than the planet, 
but decreases significantly as the satellite swarm collides 
and the dust level subsides. The planet shrinks somewhat 
as it cools, but in scattered light is nearly constant. At 
early times the planet's thermal emission peaks slightly 
short of 10 /jm, but moves to longer wavelengths over 
time. Thermal emission from the dust beyond 20 /im is 
brighter than the planet at all times. They become com- 
parable at mm wavelengths at late times, which shows 
the importance of including grain emission inefficiencies. 

This evolution therefore highlights wavelengths 
where irregular satellite clouds may be bright enough 
to be detected, both in absolute terms and relative to 
thermal and scattered emission from the planet. Unsur- 
prisingly these lie at the thermal and stellar peaks. The 
key to detection lies with rejection of starlight, the same 
issue faced by those looking for planets or circumstellar 
debris disks (e.g. Beichman et al. 2007, 2010). This re- 
jection is characterised by the star/planet contrast ratio 
that a particular instrument can detect, which is usually 



a function of angular separation between the two. While 
Figure 4 shows that the youngest dust clouds may be 
bright enough for detection by the Herschel Photodetec- 
tor Array Camera and Spectrometer (PACS) at 100 /im, 
the 7" resolution means that resolving these systems is 
difficult. We therefore focus on optical and near/mid- 
IR wavelengths. We model the detectability of satellite 
swarms using simple approximations to published con- 
trast ratios for several instruments. 

For the Hubble Space Telescope Advanced Camera 
for Surveys (HST ACS) the contrast is based on ac- 
tual roll-subtracted coronagraphic observations (Krist 
2006) and we set an absolute detection limit of 0.1 fj,Jy. 
We use predicted James Webb Space Telescope (JWST) 
Near IR Camera (NIRCAM) 4.4 /im and Mid IR Im- 
ager (MIRI) 11.4 contrast ratios from Beichman et al. 
(2010). For MIRI at 23 /im, we use the same contrast 
as at 11.4 /im for the same X/D (i.e. at twice the sepa- 
ration, Boccaletti et al. 2004). The other difference from 
the three shorter MIRI coronagraph wavelengths is that 
only a Lyot stop is offered at 23 fim so the inner working 
angle is larger. We set absolute detection limits of 68 njy, 
2.5 /iJy, and 50 ^iJy for 4.4, 11.4, and 23 /im respectively. 
These sensitivities are based on a 5 a detection in the dif- 
ference of two lhour exposures (Beichman et al. 2010). 
The contrast ratios are shown in Figure 5. 

Figure 6 shows detectability contours for the Jupiter- 
mass planet at a range of possible planetary semi-major 
axes from our A5 star at 10 pc. For each instrument there 
are two curves; one for the swarm (black curve) and one 
for the planet (grey curve). Swarms and planets in the 
space to the left of curves are detectable. The optimum 
planet semi-major axis for detection in scattered light for 
our chosen parameters (eq. 16) is drawn as a dotted line. 

Looking at the detection limits for each instru- 
ment individually, HST ACS detects scattered light from 
swarms across a wide range of the parameter space. As 
expected from equation (16), these swarms lie many tens 
of AU from the star. Detectable swarms are at larger 
semi-major axes at later times as swarms around plan- 
ets on closer orbits deplete. As expected from Figure 4, 
Jupiter-mass planets are hard to detect with ACS at any 
separation (so there is no grey HST curve on Fig. 6). 

JWST instruments are well suited to planet detec- 
tion by design so planets and swarms are detectable. The 
detection space for satellite swarms with NIRCAM cov- 
ers semi-major axes greater than 70-120 AU until a Gyr. 
Planets are detectable for a shorter time but to much 
closer separations. The difference is because the swarms 
are detected in scattered light, and planets from their 
thermal emission. At 11 fim MIRI does not detect the 
satellite swarms, but is ideal for detecting cooling plan- 
ets, whose SEDs initially peak in or near this region. The 
detection region therefore covers even more age and sep- 
aration space than at 4.4 fim. 

At 23 pm MIRI detects the cooling planet, but also 
sees the Wein side of the dust cloud for systems up to 
3 x 10 8 years old. At this wavelength the visibility space 
for MIRI lies closer than suggested by a opt because the 
dust clouds are hotter. With the exception of the vertical 
parts of curves (most notably at 23 /im), detectability 
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Figure 6. Regions in a p i vs. time space where satellite swarms 
(dark lines) and their host planets (grey dashed lines) can be 
detected with various instruments. All space to the left of each 
contour is detectable by that instrument (swarms and planets 
become fainter as they move right). The dashed line shows 
<iopt(eq. 16). The top panel shows detection contours for a 
Jupiter-mass host, and the lower panel contours for a 20 Mq 
planet. Fomalhaut b is marked by a filled circle (see §3.3). 



in this example is set by contrast rather than absolute 
sensitivity. 

Therefore, optical wavelengths appear to be the 
best place to detect extrasolar irregular satellite swarms 
around Jupiter-mass planets. Though detectable, at IR 
wavelengths the swarms are fainter than their host plan- 
ets. Discovery then involves finding whether a (presum- 
ably previously known) planet has an optical or IR excess. 

Over the parameter space in Figure 6 we find that the 
largest excess ratio (-F 8 warm/-P*) is about 2% at 100 ;im 
and 3% at 160 ^m for a young swarm at 70-80 AU. These 
excesses would not be detectable given current uncertain- 
ties with Herschel PACS photometry (Poglitsch et al. 
2010) and planet atmospheres. Given the resolution, a 
sufficiently bright swarm detected this way would be 
classed as a circumstellar disk unless it were resolved at 
shorter wavelengths or orbiting a nearby star. Such a lim- 
itation may be overcome by high resolution facilities such 
as ALMA. 

The predictions of Figure 6 for all instruments can be 
extended in several directions. For the same star moved to 
5 pc, swarms and planets at the same semi-major axis are 



more detectable because they lie at larger angular sepa- 
rations (where the instrument contrast ratio is better). 
Alternatively, a swarm around a planet at 10" separation 
around a star at 10 pc (planet at 100 AU) is about equally 
detectable at optical/near-IR wavelengths when the star 
is at 5pc (planet at 50 AU), but more easily detectable 
at mid-IR wavelengths because the dust is hotter. For 
the same range of planet semi-major axes shown in Fig- 
ure 6, the same swarms become much harder to detect 
beyond a few tens of parsecs. The difficulty arises due 
to the poorer contrast at smaller angular separations. If 
swarms are detectable beyond a few tens of parsecs, they 
must have of order an Earth-mass or more in irregular 
satellites. 

As the stellar mass decreases the detection space at 
optical/near-IR wavelengths increases due to the much 
lower stellar luminosity and consequently smaller min- 
imum grain size. However, the clouds become unde- 
tectable in the mid-IR because the Wien side of their 
spectrum is too faint. The planets themselves become 
somewhat more detectable at IR wavelengths because 
their thermal emission is the same but the star is fainter 
(and we have assumed all planetary luminosity is intrin- 
sic). 

Decreasing the planet mass can make for unde- 
tectable planets that host detectable swarms. The weak 
Hill radius dependence on planet mass means that colli- 
sional evolution is only somewhat slower (eq. 7). In ad- 
dition, the minimum grain size only depends weakly on 
planet mass (eq. 9). To highlight the relative unimpor- 
tance of planet mass on the evolution of satellite swarms, 
the lower panel of Figure 6 shows detectability contours 
for a 20 M® planet. The much fainter planet can only 
be detected at the earliest times with MIRI at 11 /im. 
Despite the 16 times decease in planet mass, the HST 
detection space is not much smaller than for the Jupiter- 
mass planet in the upper panel and still covers a wide 
range of orbits and ages. The detection space for NIR- 
CAM at 4.4 ^im is extremely small. At 23 /im, MIRI can 
still detect young swarms orbiting planets at 25 to 80 AU 
from the star. This Figure therefore shows that relatively 
low-mass planets that would otherwise be invisible can 
be detected thanks to the luminosity of their irregular 
satellite swarms. As we argue below for the specific case 
of Fomalhaut b, these swarms could be misidentified as 
more massive planets in the first instance. 



3.2.1 Discussion 

The ability of planets to capture irregular satellite 
swarms and reside at tens of AU sets the likelihood that 
any will be detected outside the Solar System. The wide 
range of planet-masses about which swarms should be 
detectable means that how capture mechanisms and effi- 
ciencies change with planet mass is important. In partic- 
ular, our prediction that swarms may be visible around 
relatively low-mass planets relies on the ability of these 
planets to capture swarms. 

Giants forming by core accretion may capture ir- 
regulars passing through their primordial envelopes (e.g. 
Pollack et al. 1979), or by "pull-down" during a phase of 
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rapid growth (e.g. Heppenheimer & Porco 1977). Giants 
that form by gravitational instability may have analogous 
capture mechanisms. Jewitt & Haghighipour (2007) sug- 
gest that the similarity of the irregular populations at 
each planet argues against gas giant specific capture pro- 
cesses for the Solar System. However, Bottke et al. (2010) 
note that because the size distributions are a result of the 
collisional evolution, they cannot be used to constrain the 
capture mechanism. 

Any planet may capture satellites during three 
and n-body interactions (e.g. Colombo & Franklin 1971; 
Agnor & Hamilton 2006). Planets may capture satellites 
as they themselves interact in the presence of a planetesi- 
mal disk. Irregulars captured earlier by gas drag and pull- 
down are supplanted by those captured during planet- 
planet interactions. Such a scenario has been proposed 
for the origin of the Solar System's irregular satellites 
within the context of the Nice model (Nesvorny et al. 
2007). The strength of this model is that it results in 
similar populations at each planet. In fact, in most model 
runs Nesvorny et al. (2007) find that Jupiter captures the 
least irregular satellites, because it undergoes the fewest 
encounters with other planets. The prospects for lower 
mass planets harbouring swarms of irregulars are there- 
fore good. 

Because we predict that swarms are most detectable 
around planets with large semi-major axes, the proba- 
bility of detection depends on the ability of planets to 
form on, or move to such orbits. The recent discovery of 
planets around HR 8799 (Marois et al. 2008) and Fomal- 
haut (Kalas et al. 2008) on such wide orbits is encourag- 
ing, not only due to their very existence, but because 
both stars also harbour planetesimal disks (Aumann 
1984; Sadakane & Nishida 1986). Planets that originate 
on closer orbits must either scatter or migrate to such 
distances. If these systems have planetesimal disks, satel- 
lite capture during the scattering and migration process 
is likely. In such a scenario, the swarm's radial extent 
(rj) cannot be much larger than half the Hill radius at 
the planetary semi-major axis where the swarm was cap- 
tured (see also §3.3.3). 

The presence or absence of irregular swarms within 
the context of discovered planetary systems should pro- 
vide information about planet formation and evolution. 
For example, the Nesvorny et al. (2007) scenario requires 
the apparently specialised circumstances of a set of giant 
planets that are destabilised and interact in the pres- 
ence of a planetesimal disk, but appears to be a robust 
way to capture irregulars. Such a scenario suggests that 
swarms may be discovered around planets in multiple sys- 
tems that harbour debris disks (such as HR8799), espe- 
cially those where scattering rather than migration has 
occurred. On the other hand, if capture of irregulars is 
associated with major planetesimal depletion events (as 
proposed by Nesvorny et al. 2007), systems with irreg- 
ular swarms may tend not to have visible debris disks. 
If primordial irregulars are formed in all planetary sys- 
tems, the properties of swarm-harbouring systems may 
be much more general. Therefore, it will be interesting 
to study the expected detection outcomes for a range of 
capture and formation mechanisms. 

As noted previously, there are effects beyond the 



scope of our model that can affect our predictions. More 
theoretical effort should be made to understand which 
are important and how they affect our conclusions, in 
particular the effect of radiation forces on dust. Efforts 
should be made to detect dust clouds around Solar Sys- 
tem giant planets, not only for their intrinsic interest, but 
to provide an empirical conversion between the mass in 
irregulars and the surface area in dust to calibrate our 
model predictions. 



3.3 Fomalhaut b 

An interesting application of our model is the planet or- 
biting the nearby star Fomalhaut, which also harbours 
a narrow ring-like circumstellar debris disk (Kalas et al. 
2008). The planet was predicted to exist based on the el- 
liptical shape of the debris ring (Kalas et al. 2005). Using 
the sharpness of the inner edge of the dust ring Quillen 
(2006) estimated the planet's orbit and a range of possi- 
ble masses between Neptune's and Saturn's. In fact the 
structure of the ring continues to provide the most strin- 
gent constraints on the planet's mass, with more recent 
modelling providing an estimate of <3 Af j up if this planet 
is the "sole sculptor" of the disk (Chiang et al. 2009). 

Chiang et al. note that Fomalhaut b need not be the 
only sculptor of the belt. This point is in part made 
because the derived orbital velocity for Fomalhaut b is 
marginally inconsistent with that expected if the disk 
and planet are apsidally aligned. With Fomalhaut b only 
observed at two epochs, this inconsistency is a minor is- 
sue at present, but highlights the possibility that another 
planet may be partly responsible for the offset and trun- 
cation of the debris disk, and Fomalhaut b may be much 
less massive. 

Another reason the mass is so poorly constrained 
at present is the lack of information about the planet's 
spectrum (Kalas et al. 2008). Thus far it has only been 
detected in two bands (0.6 and 0.8/im), the first of which 
shows a factor 3 change in brightness between measure- 
ments at two epochs separated by 2 years. Neglecting the 
potential temporal evolution, these observations and the 
non-detections in other wavebands more closely resemble 
reflected starlight than thermal emission from a planet. 

These observations lead Kalas et al. to suggest that 
the planet hosts a circumplanetary ring akin to Saturn's, 
which would extend to at least 20 Jupiter radii for an 
assumed albedo of 0.4 to recover the observed fluxes. If 
the emission is indeed scattered light from dust then the 
planet mass could be much less than 3 A/j up (see also 
Arnold & Schneider 2004). In this section we suggest that 
the observed spectrum could indeed be reflected starlight, 
scattered off the dusty debris from a circumplanetary 
swarm of irregular satellites. 

A convenient parameter space for this problem is the 
cloud semi-major axis in units of Hill radii (rj) vs. the 
planet mass (Af p i). The observed quantities constrain n, 
because the planet is not resolved and the disk size is lim- 
ited by optical depth. The planet mass will be constrained 
because a circumplanetary swarm that evolves like our 
model will not survive for the age of Fomalhaut over all 
parts of this parameter space. That is, for a swarm or- 
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biting a ~Jupiter-mass planet to be unresolved it would 
have small 77, where it would be rapidly depleted to un- 
detectable levels. 



3.3.1 Observational Constraints 

The scattered light model in Kalas et al. (2008) requires 
a total cross-sectional area of dust: 



ertot =5.8 x 10~ 4 (0.08/Q)(l/cosi) AU 2 



(24) 



where we assume Q = 0.08 (see §2.4) and cosi is the 
factor that would be required if the geometry was that of 
a flat ring inclined i to our line-of-sight, assumed to be 
~1 as the model of §2 is optically thin. 

The light from Fomalhaut b looks unresolved so the 
dust must be confined to a region smaller than the Hub- 
ble PSF FWHM of 0.53 AU. This region has an area of 
0.22AU 2 so the geometrical optical depth of the dust 
could be as low as 2.6 x 10 -3 . The resolution constraint 
means that 777-H < 0.53 so 77 < Q.6/M^( 3 , where we have 
allowed the dust cloud to be a factor of two larger in ex- 
tent to account for the concentration of brightness closer 
to the planet seen in Figure 3. As noted in §2, another 
constraint on 77 comes from the stability of circumplane- 
tary orbits, 77 < 0.5. 

Assuming that the dust is uniformly projected on 
the sky across 7r(77rH) 2 , its geometrical optical depth is 
t = o-tot/M^rH) 2 ). Since r < 1 the dust must be located 
at 77 > 0.015/Mp/ 3 , which sets a lower limit on planet 
mass of 26 x 1O _6 M0 (i.e., 6 times less massive than 
Ceres and 80 times less massive than Pluto). 

The above constraints, and the dynamical constraint 
°f -Mpi 3Mj U p from Chiang et al. (2009), are sum- 
marised by solid lines in Figure 7. The satellite cloud may 
reside anywhere in the region enclosed by the solid lines, 
which spans over seven orders of magnitude in planet 
mass, and more than two in semi-major axis. These con- 
straints are independent of our model (except perhaps 
the 0.5 Rh stability limit), and apply to any circumplan- 
etary dust population. The ~20 i?j up disk proposed by 
Kalas et al. (2008) would lie in the lower right region of 
the allowed space in Figure 7. 

Our aim is now to narrow this parameter space fur- 
ther. Because the collision rate depends strongly on 77, 
satellite clouds appear more likely to survive for the age 
of Fomalhaut around planets much less massive than 
Jupiter. 

Currently, the typically quoted age of Fomalhaut is 
200 Myr (Barrado y Navascues et al. 1997). However, re- 
cent work suggests that the age may be closer to 400 Myr 
(E. Mamajek, priv. comm.). We adopt 200 Myr as the 
age, and note below the (small) difference an older age 
makes to our model. 



3.3.2 Collisional constraints 

To find where satellite clouds survive for the age of Fo- 
malhaut, we set the collision time (1/R CC ) to 200 Myr 
(i.e. it is collision limited). It remains to substitute ap- 
propriate equations and estimates for parameters to re- 
duce Equation (7) to 77 as just a function of M p i. We 



0.010 r 



10 10 10 10 

M p , (Me.*) 

Figure 7. Constraints on the location of the proposed irreg- 
ular circumplanetary swarm (solid lines) and loci indicating 
lower limits to allowed regions with our model (dotted and 
dashed lines). The region where a satellite swarm could sur- 
vive for 200 Myr is enclosed by the dashed and dotted lines, 
and the solid resolution line (shaded grey). 



use M* = 2M , = 21L Q , p = 1500kg m~ 3 , a pi = 
115 AU, and D c = 500 km. For comparison with Table 
1, the Hill radii for Earth, Neptune, and Jupiter-mass 
planets at this distance are 0.9, 2.3, and 6.2AU respec- 
tively. We need to convert the observed <xtot into Mtot, 
for which we use equation (4). Because the effect of radi- 
ation pressure in setting D m i n varies over the parameter 
space, we also need equation (9), which gives the mass 
as Mtot = 300<7 t ot77 0,35 /Mp; 23 . The mass in satellites is 
therefore about 0.06 Mq, or 5M moon for 77 = 0.2 and 
Mpi = 10 Me, about 100-1000 times more massive than 
the initial conditions in the Bottke et al. (2010) models, 
and similar to §3.2. Such a low planet mass and higher 
stellar mass and luminosity result in D m i n ~ 300 (im, 
much higher than in the Solar System. 

Substituting these parameters into equation (7) 
yields a locus for a collision limited satellite swarm 
around Fomalhaut b that reproduces the observed ctot- 
This locus 77 = 0.27/M°! 12 is shown in Figure 7 as a dot- 
ted line. Because D c may be smaller than our assumed 
500 km, but based on the Solar System is unlikely to be 
significantly larger, the space above this line is also al- 
lowed. Thus, the collision rate in concert with the resolu- 
tion limit constrains the planet to have a maximum mass 
of about 100 M® . The minimum mass is lower than the 
mass in satellites, and therefore physically unreasonable. 

There is a further complication due to the large semi- 
major axis of Fomalhaut b, which results in low collision 



velocities t> rc i = 68M^( 3 / y^rj (in m s -1 ). A rough estimate 
of the largest object that can be destroyed by another 
of the same size (and will therefore participate in the 
collisional evolution) can be derived from equation (2) 
and X c = 1 (§2.1), giving D c « (in km). This 

estimate yields another collision rate from equation (5) 
that applies when D c is set by collision velocities 



Rc 



114 



Mtot Af* 



D° c 37 p( V a pl ) 3 M pl 



(25) 
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in years" . This equation yields another locus in the r\ vs. 
M p i parameter space where the collision rate equals the 
system age of r\ — 0.67/Mpf 46 . This limit is shown in Fig- 
ure 7 as a dashed line. Though it would appear that the 
dotted line allows swarms with D c — 500 km to survive 
for planet masses right down to 10 _3 M^, satellites this 
large cannot be destroyed at the low collision velocities 
around such low mass planets. The dashed line is a lower 
limit because D c could be smaller than the smallest ob- 
ject that can be destroyed. This limit therefore constrains 
the mass of Fomalhaut b to be more than a few Earth 
masses. 

These two loci combined with the previous con- 
straints map out a region of parameter space in Figure 7 
where a satellite swarm could survive for 200 Myr around 
Fomalhaut b. The planet is ~ 2-100 Earth masses and the 
swarm lies at 0.1-0.4 Hill radii. The swarm mass is of or- 
der a few Lunar masses (but varies with planet mass due 
to the changing D m i n , see above). The model predicts 
that a cloud of the observed luminosity could survive 
around a more massive planet, but this possibility is ex- 
cluded because Fomalhaut b is unresolved. If Fomalhaut 
is 400 Myr old, the allowed parameter space is pushed 
even closer to the resolution limit, but does not change 
Figure 7 significantly. 

Fomalhaut b is marked in the lower panel of Figure 6 
(based only on age and orbit). As expected it lies within 
the region where swarms can be detected around a 20 
planet with HST and JWST NIRCAM. 

Of course there is considerable uncertainty in both 
our model and the parameters it uses, but based on ev- 
idence that satellite swarms exist around Solar System 
planets, it at least suggests that circumplanetary swarms 
should be considered a possibility around extra solar 
planets. In the particular case of Fomalhaut b, the ques- 
tion of planet vs. cloud can resolved if new observations 
show Fomalhaut b looks like a planetary atmosphere. In 
the event that scattered light cannot be ruled out, both 
circumplanetary rings (e.g. Kalas et al. 2008) and swarms 
remain plausible options. 2 



3.3.3 Discussion 

While our model provides a possible explanation of Fo- 
malhaut b's apparently blue spectrum, the provenance 
of such a configuration is unclear. Any scenario faces the 
difficulty of explaining how Fomalhaut b came to be at 
such a large distance from Fomalhaut itself. In a core ac- 
cretion scenario, the planet presumably originates some- 
where much closer to the star (e.g. Kennedy & Kenyon 
2008), and somehow scatters or migrates to its current 
location. Because the Hill radius expands as the planet 
moves outward, the satellite cloud would need to be cap- 
tured while the planet was orbiting at least ~20 AU from 



2 Patient observers may someday find photometric phase vari- 
ations that would be caused by an optically thick circumplane- 
tary disk (Arnold & Schneider 2004). However, other methods 
will likely become available before any variation can be found 
due to Fomalhaut b's ~10 3 year period. 



the star for it to reside at rj > 0.1 now. Capture of satel- 
lites might happen as Fomalhaut b scatters off other plan- 
ets or migrates through a planetesimal disk on the way to 
its current location. The gas drag and pull-down mecha- 
nisms are unlikely to operate because our predicted mass 
of Fomalhaut b is insufficient for it to have a significant 
gaseous envelope. 

Our predicted mass for Fomalhaut b is similar to or 
less than the mass of the main debris ring itself, which 
is ~3-300 M e (e.g. Wyatt & Dent 2002; Chiang et al. 
2009). If the planet and ring masses are similar, it is un- 
likely that Fomalhaut b is responsible for truncating and 
imposing eccentricity on ring particles. Even if such a 
low-mass planet could reproduce the debris ring struc- 
ture, the relatively small chaotic zone width may require 
that the planet lie closer to the ring than observed. 

These issues do not necessarily pose a major problem 
for our model. Chiang et al. (2009) emphasise that Foma- 
lhaut b may not actually be responsible for sculpting Fo- 
malhaut's debris ring, in part because the planet's orbit is 
mildly inconsistent with the expected apsidal alignment 
of planet and ring particle orbits. However, an as yet 
undiscovered object massive enough to sculpt the debris 
ring may have more serious implications for the stability 
of Fomalhaut b. 



3.3.4 An alternative model 

The current poor constraint on Fomalhaut b's orbit al- 
lows for other interesting possibilities. For example, its 
orbit may pass through the circumstellar ring. When the 
planet is within the ring, planetesimal impacts will gen- 
erate a surrounding cloud of regolith, mantle, and plan- 
etesimal fragments. The material may free-fall back to 
the planet as in the model of Wyatt & Dent (2002), or 
in a picture more like the formation of the Earth-Moon 
system or Kuiper belt binaries some fraction of the frag- 
ments can remain in bound orbits. If Fomalhaut b spends 
a non-neglible fraction of its orbit within the ring, possi- 
ble if the planet orbit and ring are coplanar, this scenario 
provides a mechanism by which a population of objects 
bound to the planet may be built up and replenished, 
thus avoiding the requirement that the swarm must sur- 
vive for the age of Fomalhaut to be observed. This way, 
even a small amount of mass launched into orbit each 
time the planet passes through the ring will build up over 
time, eventually reaching an equilibrium state where the 
bound mass is limited by collisions. Finding this equi- 
librium mass is then simply a matter of estimating the 
rate at which mass is added to the swarm and balancing 
it with the depletion rate. Because we expect the new 
satellites to be launched into near-isotropic orbits (i.e. to 
look like an irregular satellite swarm), the depletion rate 
is calculated from the model developed in §2. 

The rate at which new satellites are added to the 
swarm is uncertain, so we adopt the approach taken by 
Wyatt & Dent (2002); an impacting planetesimal can at 
most launch its own mass of regolith and fragments to 
a significant distance from the planet. For impacts by 
large planetesimals these collisions may be more akin 
to the "graze and capture" like scenario suggested for 
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Pluto-Charon and Haumea's satellites (e.g. Canup 2005; 
Leinhardt et al. 2010). Not all planetesimals accreted by 
the planet will launch their own mass in regolith and 
fragments away from the planet, nor will all launched 
material end up in orbit. We therefore expect that the 
accretion rate must be at least a few orders of magnitude 
higher than the loss rate for the existence of a cloud to 
be feasible. 

The gravitationally focused mass accretion rate can 
be calculated from 



M a , 



Mrln 



2 TV rf ing dr r in 



1 + -fH Vrel (26) 



We assume 



I i in 



= 141 AU, dr ri n g = 25 AU, and i Iiag = 



0.026 rad (using ring parameters from Kalas et al. 2005). 
We assume the largest ring planetesimal is 500km in 
diameter, which corresponds to a ring mass of about 
75M© (Kalas et al. 2005), which is intermediate to the 
two cases considered in Wyatt & Dent (2002). The ac- 
cretion rate of a 10 planet with mass density 5000 kg 
m~ 3 is about 10~ 10 M©/year. Because the planet can- 
not spend all its time within the ring, we (arbitrarily) 
decrease the accretion rate by a factor of two. The accre- 
tion rate would be much lower for a non-coplanar orbit, 
on which the planet would spend much less time (if any) 
within the ring. For a 10 M© planet with mass density 
5000 kg/m 3 , and u rc i = 0.1 vk, Equation (26) reduces to 
2 x 10- 10 M p 4 / 3 M©/yr. 

The mass loss rate is simply the current mass times 
the collision rate, which for a swarm in collisional equi- 
librium is also the inverse of the age (equation 14). More 
generally, using Equations (2), (4), and (7), the mass loss 
rate is 



0.002 



„0.37 n l-4 



Dtt M 



1.38 /.2.27 



4.13 



2 

Otot 



(27) 



in M©/ year. For our purposes here, the mass loss rate 
is simply proportional to the square of the observed utot ■ 
That is, more mass is lost and at a faster rate for higher 
fftot or Mtot- All parameters in the large parenthesis de- 
pend on the particular model assumptions. The rate is 
independent of the maximum planetesimal size D c . This 
independence can be viewed as due to our assumption 
of a self-similar size distribution, where the mass lost is 
the same in each logarithmic size bin. If the small end of 
the size distribution is fixed, larger D c means more mass 
is available to be depleted (M to t), but collisions are less 
frequent and the largest objects are stronger. 

Evaluating Equation (27) for Fomalhaut b yields 
M loss = 5.5 x lO" 12 / (t7 3 - 4 M° 1 7 ). Using the system 
age and a planet with a collision limited swarm of 
mass 0.06 M© from the previous section yields 3 x 
10~ 1() /year, in agreement with the general expression 
for these parameters. 

Using these two expressions for the mass accretion 
and loss rates and assuming an efficiency parameter 
/acc, a swarm will replenished as fast as it decays when 
/ acc Mace = M\ OBB . Figure 8 shows where this equation is 
satisfied for different f acc . The line where the planet mass 
is too large to be consistent with the observed debris ring 
structure is omitted because a planet that passes through 
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0.001 
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Figure 8. Lines in r\ vs. M p \ parameter space where swarms 
can be replenished by accretion with different efficiency pa- 
rameters / acc (shown in legend). Though the "disk disturbed" 
limit is not drawn, we expect it to lie at a much lower planet 
mass than in Figure 7. 



ring will perturb it in a different way. This line would 
move to lower planet masses, and constrain this model. 
Another constraint on the upper planet mass would come 
from the mass above which a significant atmosphere is 
present, meaning that ring planetesimals are engulfed, 
rather than launch regolith and fragments. However, an 
atmosphere under periodic bombardment would be more 
prone to thermal escape due to heating. 

The location of the lines in Figure 8 suggest that 
replenishing a swarm from ring planetesimals is difficult. 
Even if a tenth of the accreted mass is launched into or- 
bit the planet must be more than an Earth mass, which 
seems likely to significantly perturb the main ring on a 
timescale less than 200Myr. One way to lower the planet 
mass but maintain the accretion rate is if the planet is a 
binary. The cross section for an interaction with a pass- 
ing planetesimal scales with the binary semi-major axis 
and is therefore much larger than the physical size of 
either object (Hills 1975; Leonard 1989). Such an inter- 
action does not guarantee a collision, but allows a greater 
chance of one for the same (total) planet mass. This may 
allow a binary to launch sufficient material into a swarm, 
yet also have a low enough mass to minimise dynamical 
perturbations to Fomalhaut's debris ring. 



4 SUMMARY AND CONCLUSIONS 

The purpose of this paper is to derive a simple but gen- 
eral model for collisional evolution of irregular satellites 
and apply it to planets in the Solar System and else- 
where. Our model uses the particle-in-a-box formalism 
to describe the collision rate of the largest objects in a 
satellite swarm. A model size distribution allows the mass 
evolution to be converted to the surface area in dust and 
flux densities. Though the model can be developed fur- 
ther in many ways, a comparison with the Solar System 
irregular satellites shows that our model is reasonable. 

Application of our model to the Solar System sug- 
gests that the bulk of grains may be lost to the planet 



© 0000 RAS, MNRAS 000, 000-000 



18 G. M. Kennedy & M. C. Wyatt 



or interplanetary space when the cloud is collision domi- 
nated. It may be that deposition on regular satellites only 
becomes important at later times when PR drag starts 
to dominate grain orbital evolution. Some level of dust 
must be present if the irregular satellites are still grind- 
ing down, which we suggest may be at detectable levels 
at any of the Solar System's outer planets. 

A remarkable feature of irregular satellite swarm evo- 
lution is its relative insensitivity to planet mass. Swarms 
are nearly as bright and last nearly as long around 
Neptune-mass planets as around Jupiter-mass planets. 
In contrast to extrasolar planets, we find that satellite 
swarms are most visible at optical wavelengths around 
planets that orbit at many tens of AU from their par- 
ent stars. Lower mass planets without gaseous envelopes 
cannot capture swarms by gas drag and pull-down, but 
dynamical mechanisms can still operate. There is an opti- 
mum distance for detection, which arises because swarms 
around planets on close orbits decay too rapidly and 
swarms around planets on distant orbits take too long 
to collide. 

We propose a plausible model for a satellite swarm 
around Fomalhaut b. This model provides an alterna- 
tive explanation for the spectrum of the planet being 
consistent with reflected starlight. In this picture Foma- 
lhaut b is predicted to be ~ 1-100 Earth masses, further 
illustrating that planet mass is relatively unimportant 
for the evolution of irregular satellite swarms. The al- 
lowed parameter space for the model lies very close to 
the resolution limit, so observations at higher resolution 
(FWHM < 0.53 AU or 69 mas) should test our model. 
However, the order-of-magnitude spirit of our model and 
large uncertainties in model parameters that affect our 
allowed parameter space, such as satellite size, composi- 
tion and strength, mean that the best way to test our 
model in this particular instance is to ascertain whether 
Fomalhaut b has the spectrum expected from a planet. 

We briefly outline another possibility for a swarm 
around Fomalhaut b based on speculation that it could 
pass through Fomalhaut's large debris ring. In this case, 
a swarm of objects is maintained by regolith and frag- 
ments launched by planetesimal impacts. The ~Earth 
mass planet required by such a model may be too high to 
avoid large (and unobserved) perturbations to the main 
ring. This issue may be somewhat alleviated if the planet 
is actually a binary, which would enhance the accretion 
rate and allow a lower total planet mass. 

The prospects for detecting dust created by irregular 
satellite collisions, both in the Solar System and around 
planets orbiting other stars appear good. Along with fur- 
ther characterisation of Fomalhaut b, the four popula- 
tions of irregulars on our doorstep seem a good place to 
start. 
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